Entangling flux qubits with a bipolar dynamic inductance 
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We propose a scheme to implement variable coupling between two flux qubits using the screening 
current response of a dc Superconducting QUantum Interference Device (SQUID). The coupling 
strength is adjusted by the current bias applied to the SQUID and can be varied continuously from 
positive to negative values, allowing cancellation of the direct mutual inductance between the qubits. 
We show that this variable coupling scheme permits efficient realization of universal quantum logic. 
The same SQUID can be used to determine the flux states of the qubits. 

PACS numbers: 03.67.Lx, 85.25. Cp, 85.25.Dq 



A rich variety of quantum bits (qubits) is being ex- 
plored for possible implementation in a future quantum 
computer [l| . Of these, solid state qubits are attractive 
because of their inherent scalability using well established 
microfabrication techniques. A subset of these qubits is 
superconducting, and includes devices based on charge 
0, 13 , magnetic flux 0, H, la] , and the phase difference Q 
across a Josephson junction. To implement a quantum 
algorithm, one must be able to entangle multiple qubits, 
so that an interaction term is required in the Hamilto- 
nian describing a two qubit system. For two supercon- 
ducting flux qubits, the natural interaction is between 
the magnetic fluxes. Placing the two qubits in proxim- 
ity provides a permanent coupling through their mutual 
inductance |8J|. Pulse sequences for generating entan- 
glement have been derived for several superconducting 
qubits with fixed interaction energies |fj,ll0j. However, 
entangling operations can be much more efficient if the 
interaction can be varied and, ideally, turned off during 
parts of the manipulation. A variable coupling scheme 
for charge-based superconducting qubits with a bipolar 
interaction has been suggested recently 11]. For flux 
qubits, while switchable couplings have been proposed 
previously [pj, LUl , these approaches do not enable one 
to turn off the coupling entirely and require separate cou- 
pling and flux readout devices. 

In this Letter, we propose a new coupling scheme for 
flux qubits in which the interaction is adjusted by chang- 
ing a relatively small current. For suitable device pa- 
rameters the sign of the coupling can also be changed, 
thus making it possible to null out the direct interaction 
between the flux qubits. Furthermore, the same device 
can be used both to vary the coupling and to read out 
the flux states of the qubits. We show explicitly how 
this variable qubit coupling can be combined with mi- 
crowave pulses to perform the quantum Controllcd-NOT 
(CNOT) logic gate. Using microwave pulses also for ar- 
bitrary single-qubit operations, this scheme provides all 
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FIG. 1: (a) SQUID-based coupling scheme. The admittance 
Y represents the SQUID bias circuitry, (b) Response of 
SQUID circulating current J to applied flux $ 3 for /3l = 
0.092 and I b /I c (OA5$o) = 0,0.4,0.6,0.85 (top to bottom). 
Lower right inset shows J(& s ) for same values of 7), near 
$ s = 0.45$o- Upper left inset shows 7 C versus $ s . 



the necessary ingredients for implementation of scalable 
univeral quantum logic. 

The coupling is mediated by the circulating current J 
in a dc Superconducting QUantum Interference Device 
(SQUID), in the zero voltage state, which is coupled to 
each of two qubits through an identical mutual induc- 
tance M qs [Fig. Ola)]. A variation in the flux applied to 
the SQUID, $ s , changes J [Fig. Ufb)]. The response is 
governed by the screening parameter (3l = 2LIo/<&o and 
the bias current lb, where h < I c (^s), the critical cur- 
rent for which the SQUID switches out of the zero voltage 
state at T = in the absence of quantum tunneling. In 
flux qubit experiments |l7|. the flux state is determined 
by a dc SQUID to which fast pulses of lb are applied 
to measure I C (Q S ,T). Thus, existing technology allows 
lb to be varied rapidly, and a single dc SQUID can be 
used both to measure the two qubits and to couple them 
together controllably. 

The flux qubit consists of a superconducting loop inter- 
rupted by three Josephson tunnel junctions jj,ll2]- With 
a flux bias near the degeneracy point, $o/2, a screening 



current I q can flow in either direction around the qubit 
loop. Given the tunnel coupling energy S between the dif- 
ferent directions of I q , the ground and first excited states 
of the qubit correspond to symmetric and antisymmetric 
superpositions of these two current states. Thus, the dy- 
namics of qubit i can be approximated by the two-state 
Hamiltonian 



«i = -(e?/2y«-05 4 /2)aW. 



(1) 



The energy biases e° are determined by the flux bias of 
each qubit relative to <t> /2- The tunnel frequencies 8i/h 
are fixed by the device parameters, and are typically a few 
GHz. For two flux qubits, arranged so that a flux change 
in one qubit alters the flux in the other, the coupled-qubit 
Hamiltonian describing the dynamics in the complex 4- 
dimensional Hilbert space becomes 

H = Hi® I (2) + I {1) U 2 - {K/2)a^ ® a^ , (2) 

where v 1 ' is the identity matrix for qubit i and K char- 
acterizes the coupling energy. For K < 0, the mini- 
mum energy configuration corresponds to anti-parallel 
fluxes. For two flux qubits coupled through a mu- 
tual inductance M qq , the interaction energy is fixed at 
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For the configuration of Fig. Qfa), in addition to the 
direct coupling, Kq, the qubits interact by changing the 
current J in the SQUID. The response of J to a flux 
change depends strongly on If, [Fig. Cfb)]. When I q ' 

(2) 

switches direction, the flux coupled to the SQUID, A<I>s , 
induces a change AJ in the circulating current in the 
SQUID, and alters the flux coupled from the SQUID to 
qubit 1. The corresponding coupling is 
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Re(dJ/d$ a ) Ih .(3) 



The transfer function, (dJ/d$ s ) I , is related to the dy- 
namic impedance, Z, of the SQUID via |l5j 



dJ/d$ s = iuj/Z = l/C + iuj/TZ, 



(4) 



where 1Z is the dynamic resistance, determined by Y 
which dominates any loss in the Josephson junctions, and 
C is the dynamic inductance which, in general, differs 
from the geometrical inductance of the SQUID, L. 

We evaluate {dJ /d$ 8 ) I by current conservation, ne- 
glecting currents flowing through the junction resis- 
tances: 

I b = I y + 2J cos A 7 sin 7 - 2C($o/2tt)7, (5) 



J = 7 cos 7 sin A7 - C(<I>o/27r)A7. 



(6) 



Here, I y is the current flowing through the admittance 
Y(u>) [Fig. da)], and Iq and C are the critical current and 
capacitance of each SQUID junction. The phase variables 
are related to the phases across each junction, 71 and 72, 




FIG. 2: (a) Variation of K with h for $ s = 0.45$o and device 
parameters described in text, (b) Highest achievable value of 
K 3 versus /3l evaluated at If, = 0.85/ c (0.45<l?o); ^0 (and hence 
Pi) is varied for L = 200 pH. 



as A7 = (7! — 72V2 and 7 = (71 + 72)/2- The phases 
are constrained by dAj — (7r/$o)(^s ~ LdJ). 

The expression for K s in terms of Re(9J/9$ s ) 7 
(Eq. |2J) requires the qubit frequencies to be much lower 
than the characteristic frequencies of the SQUID. This 
condition is satisfied by our choice of device parameters, 
and also ensures that the SQUID stays in its ground state 
during qubit entangling operations. Furthermore, it is 
a reasonable approximation to take the u> = limit of 
Re (dJ/d$> s ) Ib to calculate K s , so that we can solve Eqs. 
|J5J and I© numerically to obtain the working point; for 
the moment we assume Y(Q) — 0. For the small de- 
viations determining K s , we linearize Eqs. (J5J and © 
and solve for the real part of the transfer function in the 
low- frequency limit: 
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1 



I - tan 2 A7 tan 2 7 



2Li 1 



2^-(l - tan 2 A7 tan 2 7) ' 



(7) 



Here, we have introduced the Josephson inductance for 
one junction, Lj = $o/27r/ocos A7COS7. For j3l ~^> 1, 
Eq. approaches 1/L, while for (3l <C 1, 

Re(dJ/d<5> s ) h = (l/2i J )(l - tan 2 A 7 tan 2 7 ). (8) 

We see that Re(dJ/d& s )j becomes negative for suffi- 
ciently high values of lb and $ s , which increase 7 and 
A7. 

We choose the experimentally-accessible SQUID pa- 
rameters L = 200 pH, C = 5 fF, and I = 0.48 
fiA, for which j3l = 0.092. The qubits are character- 



ized by /, 



(i) 



r(2) 



0.46 11A, M qs = 33 pH, and 



M qq = 0.25 pH, yielding K /h = -0.16 GHz. Choosing 
$ s = 0.45$o, Eqs. © and Q) result in a net coupling 
strength K/h = (K Q + K s )/h that is -0.3 GHz when 
h = 0, and zero when J 6 // c (0.45$ ) = 0.57 [Fig. Efa)]. 
The change in sign of K s does not occur for all (3l ■ Fig- 
ure Hfb) shows the highest achievable value of K s versus 
Pl- We have adopted the optimal design at /?/, = 0.092. 
We also need to consider crosstalk between the cou- 
pling and single-qubit terms in the Hamiltonian. When 
the coupling is switched, in addition to dJ/d$ s being 
altered, J also changes, thus shifting the flux biases of 
the qubits. The calculated change in J as the coupler is 
switched from If, = to I(,// c (0.45$o) = 0.57 produces a 



change in the flux in each qubit corresponding to an en- 
ergy shift Sei/h = 8tijh = 1.64 GHz. In addition, when 
the qubits are driven by microwaves to produce single- 
qubit rotations, the microwave flux may also couple to 
<& s . As a result, K is weakly modulated when the cou- 
pling would nominally be turned off. A typical microwave 
drive Ci{t)/h of amplitude 1 GHz results in a variation of 
about ±14 MHz about K = 0. 

When the bias current is increased to switch off the 
coupling, the SQUID symmetry is broken and the qubits 
are coupled to the noise generated by the admittance 
Y . We estimate the decoherence due to this process by 
calculating the environmental spectral density J(y>) in 
the spin-boson model |lfj. We obtain J (to) from the 
classical equation of motion for the qubit flux with the 
dissipation from Y coupled to either qubit through J: 

J{u) = (l 2 q M*Jh) Im (dJ/d<S> s ) h . (9) 

To calculate i7(w), we linearize Eqs. © and © around 
the equilibrium point to obtain 

2 tan 7 tan A7 1 



dj 



L, 



2/Lj - 2u 2 C + iujY 



dA 1 . (10) 



For the case Y^ 1 = R. following the path to the static 
transfer function Eq. (J7J and taking the imaginary part 
in the 1ow-/?l limit, we obtain Im(dJ/d& s )i b = —uj/IZ — 
(w/4i?)tan 2 A7tan 2 7. Thus J(w) = aui, where a = 
(M^Jl/AhR) tan 2 A 7 tan 2 7, and a(I b = 0) = 0. As 
lb is increased to change the coupling strength, a in- 
creases monotonically. For the parameters described 
above and for R — 2.4 kO, when the net coupling is zero 
[I 6 /I c (0.45$o) = 0.57, Fig. Eta)] we find a = 8 x 1CT 5 , 
corresponding to a qubit dephasing time of about 500 
ns, one order of magnitude larger than values currently 
measured in flux qubits |17| . 

We now show that this configuration implements uni- 
versal quantum logic efficiently. Any n-qubit quantum 
operation can be decomposed into combinations of two- 
qubit entangling gates, for example, CNOT, and single- 
qubit gates [lg. Single-qubit gates generate local uni- 
tary transformations in the complex 2-dimensional sub- 
space for the corresponding individual qubit, while the 
two-qubit gates correspond to unitary transformations in 
the 4-dimcnsional Hilbert space. Two-qubit gates which 
cannot be decomposed into a product of single-qubit 
gates are said to be nonlocal, and may lead to entangle- 
ment between the two qubits |l4| . Since we can adjust 
the qubit coupling K to zero, we can readily implement 
single-qubit gates with microwave pulses as described be- 
low. 

To implement the nonlocal two-qubit CNOT gate, we 
use the concept of local equivalence: the two-qubit gates 
U\ and U2 are locally equivalent if XJ\ = fci£/2&2, where k\ 
and &2 are local two-qubit gates which are combinations 
of single-qubit gates applied simultaneously. These uni- 
tary transformations on the two single-qubit subspaces 



transform the gate U2 into U\. The local gate which 
precedes Us, fe, is given by fei ® ^22, where ^21(22) is a 
single-qubit gate for qubit 1(2), while the local gate which 
follows U-2, fci, is fci i0 fci2, where A: 11 ( 12 ) is a single-qubit 
gate for qubit 1(2) [l9|. Our strategy is to find efficient 
implementation of a nonlocal quantum gate U% that dif- 
fers only by local gates, k\ and &2> from CNOT, using the 
methods in [lj], and then to add those local operations 
required to achieve a CNOT gate in the computational 
basis, in which the SQUID measures the projection of 
each qubit state vector onto the z-axis. 

The local equivalence classes of two-qubit operations 
have been shown [14j to be in one-to-one correspondence 
with points in a tetrahedron, the Weyl chamber. In 
this geometric representation, any two-qubit operation 
is associated with the point [c±, C2, C3], where CNOT cor- 
responds to [7r/2,0, 0]. Furthermore, the nonlocal two- 
qubit gates generated by a Hamiltonian acting for time 
t can be mapped to a trajectory in this space [lj]. If 
K is increased instantaneously to a constant value, the 
trajectory generated by Eq. (EJ is well described by the 
following periodic curve 



[01,02,03] = [Kvt/h, p \sintot\ ,p jsinwij 



(11) 



Here, p is a function of the system parameters, v = 
e1eyAE 1 AE 2 , and u = (A£i - AE 2 )/2h, where AE, = 
[(e°) 2 + Sf] 1 / 2 is the single-qubit energy level splitting. 
Independently of p, this trajectory reaches [ir/2, 0, 0] in a 
time tx = nix J 10 when the coupling strength is tuned to 
K = hcu/2nv, with n a nonzero integer. 

While this analytic solution contains the essential 
physics, it is an approximation and does not include vital 
experimental features, in particular, crosstalk and the fi- 
nite rise time of the bias current pulse. To improve the 
accuracy, we perform a numerical optimization using Eq. 
(|llfl as a starting point, then add these corrections. We 
use tunnel frequencies S±/h = 5 GHz and 82 /h = 3 GHz, 
and include the shifts of the single-qubit energy biases 
due to the crosstalk with K s in Eq. (jl 1|) by adding a 
shift Set proportional to K. We account for the rise and 
fall times of the current pulse by using pulse edges with 
90% widths of 0.5ns [see K(t) in Fig. EJ. We numeri- 
cally optimize the variable parameters to minimize the 
Euclidean distance between the actual achieved gate and 
the desired Weyl chamber target CNOT gate. We find 
K/h = -0.30 GHz, e\/h = 8.06 GHz, e° 2 /h = 2.03 GHz, 
and tx = 8.74 ns; tx is the time during which the qubit 
coupling is turned on. 

As outlined above, to achieve a true CNOT gate we still 
have to determine the pulse sequences which implement 
the requisite local gates that take this Weyl chamber tar- 
get U2 to CNOT in the computational basis. Local gates 
may be implemented by applying microwave radiation, 
6j(t), which couples to aj , and is at or near resonance 
with the single-qubit energy level splitting AEi . We note 
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cancellation of the direct mutual inductive coupling Kq 
between the qubits so that the net coupling K can be 



5 10 15 20 25 30 

FIG. 3: Pulse sequence for implementing CNOT gate. En- 
ergy scales in GHz. Total single-qubit energy bias ei(t) = e° + 
t-iif) + 5ei(t), where microwave pulses £1,2 (£) produce single- 
qubit rotations in the decoupled configuration; crosstalk mod- 
ulation of K(t) is shown (see text). The bias current is pulsed 
to turn on the interaction in the central region. 



that the single-qubit Hamiltonian of Eq. JQ) driven by a 
resonant oscillating microwave field does not permit one 
to use standard NMR pulses, since the static and oscil- 
lating fields are not perpendicular, but rather are canted 
by an angle tan _1 ((5j/e"). To simplify the pulse sequence, 
we keep ej 2 constant at the values used for the non-local 
gate generation. This imposes an additional constraint 
on the local gates: to generate a local two-qubit gate 
fci = fcn £g> fci2, the two single-qubit gates ku and k\2 
must be simultaneous and of equal duration. We satisfy 
this constraint by making the microwave pulse addressing 
one qubit resonant and that addressing the other slightly 
off-resonance. Using this offset and the relative ampli- 
tude and phase of the two microwave pulses as variables, 
we can achieve two different single-qubit gates simulta- 
neously, leading to our required local two-qubit gate. 

The resulting pulse sequences for K and e\^ are shown 
in Fig. |3 The gate has a maximum deviation from CNOT 
in the computational basis of 1.6% in any matrix element. 
This error arises predominantly from the cross-coupling 
of the microwave signals for the two qubits and the weak 
modulation of the K = state of the coupler during 
the single-qubit microwave manipulations. While small, 
this error could be reduced further by performing the nu- 
merical optimization with higher precision or by coupling 
the microwave flux selectively to each of the qubits and 
not to the SQUID. The total elapsed time of 29.35 ns is 
comparable to measured dephasing times in a single flux 
qubit H2I. 

In summary, we have shown that the inverse dynamic 
inductance of a dc SQUID with low (3l in the zero- voltage 
state can be varied by pulsing the bias current. This 
technique provides a variable-strength interaction K s be- 
tween flux qubits coupled to the SQUID, and enables 
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FIG. 4: Chain of flux qubits with intervening dc SQUIDs 
arranged to provide both variable nearest neighbor coupling 
and qubit readout. 



switched from a substantial value to zero. By steering a 
nonlocal gate trajectory and combining it with local gates 
composed of simultaneous single-qubit rotations driven 
by resonant and off-resonant microwave pulses, we have 
shown that a simple pulse sequence containing a single 
switching of the flux coupling for fixed static flux biases 
results in a CNOT gate and full entanglement of two flux 
qubits on a timescale comparable to measured decoher- 
ence times for flux qubits. Furthermore, the same SQUID 
can be used to determine the flux state of the qubits. This 
approach should be readily scalable to larger numbers of 
qubits, as, for example, in Fig. 0J 
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